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Abstract 



Taking into account the terms ~ m 2 in the effective chiral lagrangian , we show 
that, at 9 = ir, the theory with 2 light quarks of equal mass involves two degenerate 
vacuum states separated by a barrier. For Nf = 3, the energy barrier between two 
vacua appears already in the leading order in mass. This corresponds to the first 
order phase transition at 9 = it. The surface energy density of the domain wall 
| separating two different vacua is calculated. In the immediate vicinity of the point 

9 = it, two minima of the potential still exist, but one of them becomes metastable. 
. The probability of the false vacuum decay is estimated. 

' : 1 Introduction 

Ph. 

<D . 

It is very well known that in QCD with Nf massless quarks, chiral symmetry SU(Nf) <8> 
>• . SU(Nf) is spontaneously broken down to SUy(Nf). (This is an experimental fact for 
Nf = 2, 3. Probably, no spontaneous chiral symmetry breaking occurs for large enough 
number of flavors Nf ~ 8 — 10 |JTJ] which will be of no concern for us here). Spontaneous 
symmetry breaking means that the order parameter < q^c/'l > (i,j are flavor indices) can 
acquire an arbitrary direction in the flavor space. Massless Goldstone particles appear. 
If the free quark masses are not zero, the axial chiral symmetry is broken explicitly and 
the minimum of the energy functional corresponds to a particular flavor orientation of 
the condensate. Goldstones acquire small masses. In the real World with m u ~ 4 MeV, 
m d w 7 MeV, m s ~ 150 MeV, and 9 = 0, the vacuum state is unique. 

It is interesting to study also other variants of the theory with different values of 
masses and 9. It was known for a long time @, [| that in the theory with equal light 
quark masses and 9 = tt, there are two degenerate vacuum states. This is best seen in 
the framework of the effective chiral lagrangian describing only the light pseudogoldstone 
degrees of freedom. In the leading order in mass, the effective potential is 



V = -S Re 



1 Pcrmanent address: ITEP, B. Cheremushkinskaya 25, Moscow 117218, Russia. 
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where U = exp{2i(p a t a / F n } (0 a are pseudogoldstone fields), M. is the quark mass matrix 
and S is the absolute value of the quark condensate. 

Suppose Nf — 3, A4 — ml, and 9 = 0. The minimum of the energy is achieved at 
U = 1. For 9 = 7r, there are two different minima with U = 1 and [/ = e 27ri ' 3 l. They 
are separated by the energy barrier. The appearance of two vacuum states corresponds 
to spontaneous breaking of the CP-symmetry by the Dashen mechanism 0. 

The situation is, however, more confusing for Nf = 2. The trace of a £77(2) matrix is 
always real which means that, at M. = ml and 9 = tt, the potential ([I]) does not depend 
on U at all. That would mean that the explicit breaking of the chiral symmetry is absent, 
pions are massless at this point and the phase transition would be of the second rather 
than of the first order. We understand, however, that the chiral symmetry of the original 
QCD lagrangian is broken explicitly by the quark mass term for all values of 9 including 
9 = tt, and the situation looks paradoxical. 

This question was left undiscussed in the original papers [0, [| , and only comparatively 
recently it was shown how the paradox is resolved ||. To this end, one should to take 
into account the terms of order ~ m 2 in the effective chiral lagrangian. If doing so, the 
continuous vacuum degeneracy at 9 = tt is lifted, and we obtain again only two vacuum 
states separated by a barrier. 

In what follows, we confirm this finding. The main aim of the paper is to bring the 
analysis of Ref. || into contact with standard chiral theory notations and wisdom and to 
perform some quantitative estimates both for Nf = 3 and Nf = 2 for the height of the 
energy barrier, the surface energy density of the domain walls interpolating between two 
degenerate vacua at 9 = tt, and for the decay rate of metastable vacuum states at the 
vicinity of the phase transition point. 



2 Three flavors 



Let us discuss first in some details the case Nf = 3 where no complications due to higher- 
order terms arise. By a conjugation U — ► VUV', any unitary matrix U can be brought 
into the diagonal form U = diag(e m , e l/3 , e~ i(a+/3 )). When Ai = ml, a conjugation does 
not change the potential ([!]). For diagonal U, the latter acquires the form 



U(a,/3) 



-mS 




cos (« + /?+- 



The function U has six stationary points: 



a 



(3 = 0, II 



a 



P 3 ' 



III : 



a 



(2) 



IV : a = (3 = -- 
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7T, IVa : a 



29 29 
-a -(3 = h7r, IV6: (3 = -a- (3 = -- 



71 



3 ' '3 

2 One could, of course, anticipate it. For 8 = and in the leading order in quark mass, the mass 
of pions is given by the Gell-Mann - Oakes - Renner relation F%M% = (m u + md)S. A theory with 
m u = md and 9 = tt is equivalent to the theory with m u = —nid and 9 = (only 6* p hy S = 9 + arg(detA'J) 
is relevant). In that case, the pions seem to stay massless. However, the Gell-Mann - Oakes - Renner 
relation is true only in the leading order in m q . Higher order corrections bring about a nonzero mass to 
pions. 



(3) 
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The points IVa and IVb are obtained from from IV by Weyl permutations and their 
physical properties are the same. Actually, we have here not 3 distinct stationary points, 
but the whole 4-dimensional manifold SU(3)/[SU(2) ®C/(l)] of the physically equivalent 
stationary points related to each other by conjugation. The values of the potential at the 
stationary points are 

n + 2n „ 9-2tt 

Ei = — 3mS cos -, En = — 3mS cos , Em = — 3mS cos , 

3 3 3 

£iv = m£ cos 9 (4) 

Studying the expressions (f|) and the matrix of the second derivatives at 8 = it, one can 
readily see that i) the points I and III are the degenerate minima; ii) the point II is 
the maximum, and Hi) the points IV are saddle points. When 9 is slighly less than tt, 
I is a global minimum while III is still a minimum, but of a local variety. The latter 
coalesces with all three saddle points at 9 = tt/2. At this point the eigenvalues of the 
second derivative matrix pass zero and, at still lower values of 9, a metastable minimum 
does not exist. When we instead make 9 larger than tt, the picture is symmetric, only the 
minima I and III change their roles. At 9 = 0, the picture is exactly reversed compared 
to what we had at 9 = tt: there is one global minimum I, two degenerate maxima II and 
III, and a surface of saddle points IV. 

Fig. 1 illustrates how the stationary points of the potential are moved when the 
vacuum angle is changed. One can show that metastable vacua are absent at 9 = also 
with physical values of masses. 

Two degenerate vacua at 9 = tt are separated by the domain wall. To find the profile 
of this wall, we have to restore the kinetic term ^-Tr^^Ud^W} in the effective lagrangian 
and seek for the field configurations depending only on one spatial coordinate x with the 
boundary conditions U{— oo) = 1, U(oo) = e 2t7T ^ 3 and realizing the minimum of the 
energy functional 

E = Ao = A J dx f -Z-Tr{d x Ud x rf} - mSRe [e te/3 Tr J (5) 

(A is the total area factor). In our case, it suffices to seek for the solutions in the class 
U = diag(e ia ^, e ta<yX \ e -2ia W). Introducing 7 = a — tt/3, subtracting the vacuum energy 
and using the Gell-Mann - Oakes - Renner relation, the expression (El) is rewritten as 



a = 3Fi / dx 



1 12 M 7 

-7 H I cos 7 

2' 3 V 2 



(6) 



The corresponding equations of motion with the boundary conditions 7(±oo) = ±7r/3 
can be readily integrated. The first integral is 

7' = M 7r ^273(cos 7 -i) 

Integrating it further, we obtain the solution 

E 2 + AE + l 

COS7 = 2(£P + E + 1) (7) 



3 



E 




Figure 1: Stationary point of E(a,/3) for different 9. Solid lines are the minima, dotted 
lines are the maxima, and dashed line are the saddle points 
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where E = exp{M n x/\/2}. The solution ([/]) is centered at x = where it passes through 
the saddle point IV. There are, of course, other solutions obtained from Eq.flT]) by a shift 
of x. Also we could have chosen the Ansatze U = diag(e ta ^ x \ e~ 2ia ( x \ e m ^) or U — 
diag(e _2m ^\ e ta ^ x \ e ta ^) and obtain two other wall solutions (with the same properties) 
passing through the saddle points IVa and IVb. The wall surface tension is 



a 



9 „ r°° EdE r ( 7T \ „ 

FlM n / = 3V2 1 - — = M n F 2 



V2 n ""Jo (E 2 + E + l) 2 \ 3^3 



Suppose now that 9 = tt + <p with < \<p\ <C 1. The energies of the vacua are not 
degenerate anymore but are splitted apart by the value 

AS w mSv^l^l (9) 

A metastable vacuum should decay with the formation of bubbles of the stable phase. 
The quasiclassical formula for the decay rate per unit time per unit volume was derived 
in §: 

T oc exp{-^7r 2 -^— 1 (10) 



2 (AS) 



where a is the surface tension the bubble. Substituting here Eqs.([| |^) [ Strictly speaking, 
at 6 7^ tt, the bubble surface tension does not coincide with Eq.(Q) but is somewhat less 
going to zero at 6 = n/2 or 9 = 3n/2. But at small \<f>\, the expression @ is correct], we 
obtain f\ 

( CF 2 1 



with 

C = 2 4 -3 5 v / 3tt 



TT 



3^3. 

Note that the numerical factor C in the exponent is tremendously large, i.e. lifetime 
of metastable states would be tremendously large (much larger than the lifetime of the 
Universe) for almost all 9 in the interval (|, not too close to its boundaries (where 
metastable states disappear). It is a real pity that such a beautiful possibility is not 
realized in Nature 0. 

3 Note the difference with the rough estimate lnT ~ — F% / (M%\(f>\ 2 ) for the same quantity in the 
Witten's paper gj. First, one has \4>\ 3 rather than \(j)\ 2 in the denominator and, second, a huge numerical 
factor pops up. 

4 Recently, Halperin and Zhitnitsky argued the existence of metastable states in the real QCD at 
= 00]. However, their arguments were based on a particular model form of the effective potential 
incorporating also glueball degrees of freedom and involving certain cusps. The status of this potential 
is not quite clear by now. 
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a) 



b) 



Figure 2: Connected and disconnected contributions to the quark current correlators. 

3 Two flavors 

As was mentioned before, we have to take into account here the terms of higher order 
in mass in the effective potential. For Nj = 2, it is convenient to make benefit of the 
fact that SU(2) £g> SU(2) = 0(4) and to use the 4-vector notations so that U = U^cr^ = 
Uq + iUi<Ji, U 2 — 1. The 2x2 complex mass matrix involves 8 real parameters which is 
convenient to "organize" in two different isotopic 4- vectors Q 

X, = ^ (Re Tr{W e/2 }, Im Tr{Me w / 2 a t }) 

TT 

Xf, = ^ (im Tr{Me^ 2 }, -Re Tr{Me»' 2 <Ti}) (12) 

7T 

In the second order in x, %•> the most general form of the potential is M 

v(u,) = -f 2 (xM - i 3 (xM 2 - i 7 (xM 2 (is) 

where l^j are some dimensionless coefficients (the coefficients £1,2,4,5,6 multiply the struc- 
tures involving the derivatives of the field U in the effective lagrangian). The term 
{XfiVfi){XfiUn) is not allowed because it would lead to CP breaking even at 6 = 0. 
The term ~ (x^U^) 2 can always be neglected compared to the leading one for small 
masses and is not interesting. On the contrary, the term involving I7 has a different 9 - 
dependence and, for M. = ml and 9 ~ tt when the leading term vanishes, determines the 
whole dynamics. 

Before proceeding further, let us try to extract an information on the numerical value 
of the constant I7. The following relation belonging to the same class as the well-known 

5 We use the notations of Ref. || . 
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Weinberg sum rule and derived in Ref. H is very useful: 




d x 



< S l (x)S k (0) > -5 lk < P u (x)P u (0) > 



(14) 



where S l = qa l q and P° = iqy'q. Let us calculate the right-hand side of Eq. ( |i~4"D as an 
Euclidean functional integral. Let us first calculate the correlators in a particular gauge 
field background. Only the connected quark diagram depicted in Fig. 2a contributes to 
the correlator of the scalar isovector densities. The pseudoscalar correlator receives also 
a contribution from the disconnected graph in Fig. 2b. The solid lines in Fig. 2 stand for 
the quark Green's functions in a particular gauge field background for which we use the 
spectral decomposition 



G A (x,y) = Yl 



m — iX n 



(15) 



where X n are the eigenvalues of the massless Euclidean Dirac operator in an external 
gauge field and ip n {x) are its eigenfunctions. All non-zero eigenvalues are paired: for any 
eigenfunction ip n (x) with non-zero eigenvalue X n , i> n (x) = j 5 ip n (x) is also an eigenfunction 
with A n = — A n . There are also zero modes. For each flavor, their number coincides with 
the topological charge v = (1/327T 2 ) / d^xG^G^ of the gauge field configuration. 

Now, for large Euclidean volumes mVE ^> 1, zero mode contribution is irrelevant 
for the connected graphs (cf. the discussion in Ref. 0). On the other hand, the only 
contribution in the disconnected graph for the pseudoscalar correlator is due to zero 
modes. This contribution is very large oc v 2 /m 2 and is of a paramount importance. 

pairing together 



Plugging in the Green's functions (ITq ) in the correlators in Eq.(|l^ 
positive and negative A, and integrating over gauge fields, we obtain for large volumes 



In 



1 

V 



E 



71V 



-A? 



m 2 + X 2 ) 2 



+ E 



rrr 



+ K, 



< v A > 



(16) 



where J2' n means the summation over positive eigenvalues only and the symbol < . . . > 
stands for the gauge field averaging. The first term in the R.H.S. comes from the scalar 
isovector correlator, the second term is the contribution of the connected graph to the 
pseudoscalar correlator, and the last term is due to the disconnected graph. 
Introducing the spectral density 



P(A) 



(17) 



the relation ( jig) is rewritten as 

Ft 
T 2 



7 



m 



p(X)dX 
o (A 2 + m 2 f 



< v z > 
2m 2 V 



This relation belongs to the same class as the famous Banks and Casher relation for the 
fermion condensate 



2m 



p{X)dX 
A 2 + m 2 



vrp(O) + 0(m) 



(19) 
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For some more examples, see Refs. [Hi [LI 



Both terms in (fig) are singular oc 1/rra in the chiral limit. The singularity cancels out, 
however. Indeed, the leading infrared contribution in the first term in square brackets is 
7rp(0)/(4m) ~ E/(4m) which is the same as for the second term due to the known result 
for the topological succeptibility in a theory with light quarks of the same mass (T2| 



\<v 2 >= ^ + 0(ra 2 ) (20) 

The absense of singularity in the isoscalar pseudoscalar correlator means that the corre- 
sponding meson is massive: U(l) problem is resolved by the 't Hooft mechanism due to 
fermion zero modes in topologically non-trivial gauge backgrounds, l-j is given thereby 
by a constant term ~ 0(1) which is left out after the cancellation of singular terms. This 
constant is completely determined by the term oc m 2 in < v 2 >. Indeed, for Nf = 2, the 
spectral density is analytic at A = 0: p(A) = p(0) + /iA 2 + . . . [|10|. The extra infrared 
contribution is of order 0(m) and can be neglected. R. We finally obtain 



lj = lim 



o4S 2 



2 < v 2 > /(mV) 



m 



(21) 



What can be said about m-dependence of the topological succeptibility in the next-to- 
leading order in mass ? Let us first see what happens in a theory with two light quarks 
embedded in the theory involving also the third quark which is much more massive but 
still light enough (like it is the case in the real World). The topological succeptibility in 
the theory with 3 light quarks m u = m d = m <C m s < /ihadr is given by the expression 



< v 2 > mm.S mE m 2 S o. . . 

V m + 2m s 2 4m s 

i.e. the second term is negative [] which means that l>j as given by Eq.(^Tj) is positive. We 
obtain 

h = = *m w 5 ' 10-3 (23) 



v 



The same result could be obtained in a more direct way if saturating the pseudoscalar 
correlator in Eq.(|14]) by the t] - meson pole 0. 

The estimate (|2~3D is quite good for the real QCD. What we are interested in here, 
however, is a hypothetical theory with 9 ~ it and just two light flavors. Remarkably, 
an analytic result for lj can be obtained also in this case if the number of colors iV c is 
assumed to be large. For large N c , the axial U(l) symmetry is almost not affected by 

6 Note also that the spectral integral in Eq.([l8]) involves a logarithmic ultraviolet singularity at large 
A where the spectral density is the same as for free fermions p(A) ~ A 3 . It is multiplied, however, by m 2 
and can be dropped out by that reason. 

7 That is quite natural, of course. The presence of an extra light quark brings about the suppression 
of large v due to the extra factor m v s in the fermion determinant. 
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the anomaly which means that rj' - meson is relatively light: M3 ~ (J^ aJ i T /N c . We can 
saturate now the pseudoscalar correlator by the rf pole to obtain 

F 2 

h = —^r (24) 
7 2M 2 1 ; 

The same result can be obtained via the relation fl2"l|) . For large N c and Nf = 2, the 



topological succeptibility is known to be ITS 



1 9 mrS mS m 2 X 2 ^. o N . . 

v <v >= 27T^ = T " — + 0(m) (25) 

where r is the topoligical succeptibility in the pure Yang-Mills theory. Again, the term 
~ m 2 in < v 2 > /V is negative which leads to the positive 1? which coincides with ( f2"4"|) 
due to the relation 

F 2 M 2 = At (26) 

which holds in the limit m —>■ 0. We assume that l 7 is positive also for small number of 
colors down to iV c = 2. Indeed, in the limit m — > oo, the topological succeptibility should 
coincide with r. It is natural that the series in m for small masses [the analog of Eq.(E5|)] 
should have alternating signs. We cannot, unfortunately, formulate this statement (the 
positiveness of l?) as an exact theorem though our suspicion is that such a theorem can 
somehow be proven. 

We are ready now to discuss the vacua dynamics in the region 9 ~ n. Assume 
U = diag(e ia ,e" ia ). The potential (O) (with l 3 = 0) is 



9 9 ( Yj \ 2 

Via) 2mE cos - cos a — 4/ 7 m 2 sin 2 - — cos 2 a (27) 

2 2 \F*J 

Defining <fi = 9 — tt, it can be rewritten for small |0| as 

V(a) = mE(f) cos a — 4/ 7 m 2 ( — J cos 2 a (28) 



If 



F 4 

TT 

the function (|28| ) has four stationary points: 

I a = with i?i = m£</> — 4/ 7 m 2 



<0*=^, (29) 



A 2 J 



II a = ii with i?n = —rallcj) — 4l 7 m ^— J 

</>F 4 FV 2 
III, Ilia a = ± arccos „, with £ in = -^r~ (30) 

8( 7 mE I6/7 
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E 




Figure 3: Stationary point of E(a) for Nf = 2 in the region of small \<p\ = \9 — 7r| . Solid 
lines are minima and dotted lines are maxima. 

Again, the points III, Ilia are related to each other by Weyl symmetry and we have 
actually the whole surface SU{2) /U{1) = S 2 of the equivalent stationary points. Studying 
the second derivatives d 2 V/da 2 for the branches fl30D , one readily sees that, in the region 
(p9|), the points I and II present local minima (for <fi < or 9 < ir the absolute minimum 
is I while II is a metastable state; for 9 > 7T, it is the other way round). The points III 
are degenerate maxima. The picture is depicted in Fig. 3. Physically, it is exactly the 
same as in the case Nf = 3 and we have a first order phase transition. Only the width 
of the region in 9 where two local minima coexist is much more narrow than in the case 
Nf = 3 and goes to zero in the chiral limit m — > 0. 

Note that with negative l?, the picture would be reversed compared to that in Fig. 3 
so that we would have in the region (|29| ) a surface of degenerate minima at the points III. 
That corresponds to spontaneous breaking of flavor symmetry SUv{2) || (when 9 ^ 0, 
the contributions to the partition function coming from topologically nontrivial sectors 
are not positively defined, and the Vafa-Witten theorem |14| which prohibits spontaneous 
breaking of vector flavor symmetry at 9 = does not apply). Besides, at 9 7^ ix we would 
have the maxima of inequal height at the points I and II. Instead of one first order 
phase transition at 9 = 71 we would have two consequent second order phase transitions 
at 9 = 7r ± 0*. We find this picture rather unnatural and consider it as an additional 
argument why lj should be always positive. 

It is not difficult now to perform the same program as for Nf = 3 and to find the 
surface energy density of the domain wall at 9 = 71 and the decay rate of metastable 
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vacua at the vicinity of 6 = tt. The wall configuration and the energy density at 6 = it 
are obtained by minimizing the functional 



a 



^-(d x a) 2 ~2l 7 m 2 (^\ (cos 2a -1) 



dx 



(31) 



with the boundary conditions a(— oo) = 0, a(oo) = tt. The equations of motion have a 
simple solution 



0(2:) = 2arctan 
The surface energy density is 



exp { \/8l 7 —x 



mE 



32/ 7 



(32) 



(33) 



It is much lower numerically than in the case Nf = 3 and goes to zero in the chiral limit. 
The rate of metastable vacuum decay at \<f>\ <C 0* is estimated as 



r oc exp{-12V/?-^ 



(34) 



4 Schwinger model 

It is very instructive to compare the situation in QCD4 with what happens in QED2 
(the Schwinger model) with two light fermions of equal mass. The model was extensively 
analyzed in Ref.[|]j|. For zero mass, it is exactly solvable. When the mass m is not zero, 
but much less than the gauge coupling constant g (which carries the dimension of mass in 
two dimensions), a systematic expansion in the small parameter m/g ( typically, in some 
fractional powers thereof) can be built up. 

The model is exactly equivalent to the following bosonized model 

1 1 n 2 ( \ 2 

C bos = -(DM 2 + -id,,^) 2 - — ((£+- | + C cos cos v 7 ^- (35) 



57T 



where 0+ is the heavy field and </>_ describes light degrees of freedom. Assuming the 
"conformal" normalization for the Euclidean correlator 



<: cos. 



.(x) : : cos/?0_(O) :> c =o 



1 

2 



/3 2 /2vr 



+ \X 



/3 2 /2vr 



X 



and also 

<: cosf3<f) + (x) :>c=o = 0) 
the constant C was calculated to be |16| 

m ^l/2 e7 /2 2 3/4 



C 



71 



5/4 



(36) 



(37) 



(3? 



3 We have changed the convention for C compared to m& by a factor of 2. 
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where 7 is the Euler constant. When m = 0, also C = so that light and heavy degrees 
of freedom decouple and we have just free massive Schwinger boson with the mass 



2g 2 



71 



(39) 



and a sterile massless particle. When m/0, light "quasi-Goldstone" degrees of freedom 
acquire mass and start to interact with the heavy ones. One can write the effective 
lagrangian for the light degrees of freedom which has largely the same status as the 
effective chiral lagrangian in QCD. To lowest order, we can just freeze dj + = Q/\f%7~x, and 
the lagrangian reads 



\{d^f + C coJ- 



cos v 2nd)-. 



(40) 



It is nothing else 
the vacuum energ 
Not surprisingly, 

these "pions" f\ is 
When 6 ~ n, 
are in a position 
parameter m/g. 



as the Sine-Gordon model. It is exactly solved which allows one to find 
;y, the fermion condensate of the original theory, and the mass spectrum, 
the lowest states form an isotopic triplet. A characteric mass scale of 

\2/3 



(m 2 g cos 2 



1/3 



See Ref . [161 for the exact calculation. 



- (Ccos^ .. v „„ y _ 2/ 

the leading term in the effective potential vanishes, however, and we 
to take into account higher order corrections in the Born-Oppenheimer 
It is convenient to rewrite the lagrangian (1351) in terms of x = 0+ — 



C 



+ 



c 



I— 6 I—' 

cos - cos v 2tvx — sin - sin v 2nx 



: cos v 27T0 



(41) 



The correction ~ C 2 oc m 2 in the effective potential is given by the expression 

C 2 r—,, , 



AV(x) 



COS V O7T0_ (X) 



(7 f 

cos 2 - J d 2 y\y\ < cos V2n~x(y) cos V2irx(fy > 



'0=0 



+ sin 2 - J d 2 y\y\ < sin V2n~x{y) sin V2ttx{^) >c=o 



(42) 



This formula has the same meaning as Eq.(6.30) in the paper by Coleman 
used accurately the conformal fusion rules 



gWW . x 



3 W(o) 



+ . . . 



We only 



(43) 



neglected the operators of higher dimension in Eq.(^) (they bring about the corrections 
of still higher order in m), omitted an irrelevant additive constant, and took into account 

9 The only essential difference with the pions of QCD are that they are not true goldstones and decouple 
in the chiral limit. This is related to the Merman- Wigner-Coleman theorem JltJ forbidding the existence 
of massless interacting particles and, thereby, a spontaneous breaking of a continuous symmetry in two 
dimensions. 
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Figure 4: Some graphs contributing in the effective lagrangian in the second order. Crosses 
stand for the insertion of light vertex cos(\/27r0-) and solid lines describe heavy Schwinger 
bosons. The graphs with even number of the lines contribute in k and the graphs with 
odd number in k±. 



all loops of the heavy field x( x ) drawn in Fig. 4 (Coleman only took the first graph 
which is much similar to the graph with the rf exchange which saturates the pseudoscalar 
correlator and gives the leading contribution in l 7 in QCD 4 in the large N c limit. It is 
good enough for a qualitative estimate, and, as we will soon see, the account of other 
graphs only brings about a certain numerical factor which is very close to 1.). 

Expanding the integrand in x{ x ) an d disregarding the "tadpole" contributions 

involving the correlators in coinciding points [they vanish due to the convention (|37|)1, 
and substituting the free massive boson Green's function 



< X(a0*(0) > 



e tkx d 2 k 1 



k 2 + ii 2 4tt 2 2tt 



— K Q {nx) 



(44) 



we finally obtain 



1 6 

£ cff = -j[dv<P-) 2 + Ccos- : cos V27i<j)- : + 



7r 5 / 2 C 2 
V8g 3 



o 6 . n 

Kn COS h K\ Sill - 

2 2 



: cosv 7 ^^ : +0(C 3 ) 



(45) 



where 



Kq 



-- / z 2 dz{cosh[K (z)] - 1} = .163... 
Jo 

k x = / z 2 dzsmh[K (z)] = 1.604... 
Jo 



(46) 



The constant kq is an exact counterpart of the constant I3 in the effective potential (|13D 
for QCD4. The constant k% is an exact counterpart of I7, the quantity of primary interest 
for us here. We see that k% is positive (and that serves as an additional argument why Z7 
should be positive in QCD4). Note that if plugging in just K Q (z) instead of sinhfiTo^)] 
in the integral (|46"D for K\ [that corresponds to taking into account only the graph with 
one heavy particle exchange), we would obtain Ki = 7r/2 instead of 1.604. The effect of 
higher loops increases the value of k± just by ~ 3%. 
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Even the analog of the domain wall at 9 = tx exists in this nice model. In two 
dimensions, a domain wall is just a particle. On the classical level, it is a Sine-Gordon 



soliton corresponding to interpolating between the points 0_ = and 0_ = y7r/2. Soliton 
and antisoliton form an isotopic doublet. 

It would be very interesting to analyze the quantum problem and to find the mass 
spectrum and other characteristics of the model to the order ~ C 2 . Our impression is that 
it is not so easy to do at 9 = n : the problem is that (3 = \^8tt is exactly a boundary value 
of the coupling. The Sine-Gordon theory with (3 > \^8tt is sick: the hamiltonian does 
not have a ground state etc. This displays itself in the fact that the factor multiplying 
: cos \/87r0_ : in Eq.(^) is dimensionless and we do not have a mass parameter out of 
which the soliton mass could be composed. Probably, at 9 = it, the Born-Oppenheimer 
expansion breaks down in this case and one has to analyze the full lagrangian (|3^) by 
approximate methods (cf. recent |18|). The lagrangian (|5) may still be taken at the face 



value at other values of 9 and could allow to find corrections in mass to the exact results 
of Ref. |16[ at 9 = etc. But this is beyond the scope of the present paper. 
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